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IMAGING MARKERS IN MUSCULOSKEIiETAL DISEASE 

By 

Siau-Way Liew 
Costas Tsougarakis 
Claude D. Arnaud 
Philipp Lang 
Daniel Steines 
Barry J. LInder 

Technical Field 

[001] This invention relates to using imaging metliods to aid in drug discovery and 

drug deveiopment. This invention also relates to using imaging methods for diagnosis, 
prognostication, monitoring and management of disease, particularly where that disease 
affects the musculoskeletal system. This invention identifies novel imaging markers for 
use in diagnosis, prognostication, monitoring and management of disease, including 
musculoskeletal disease. 

Background 

[002] Osteoporosis and osteoarthritis are among the most common conditions to 

affect the musculoskeletal system, as well as frequent causes of locomotor pain and 
disability. Osteoporosis can occur in both human and animal subjects (e.g. horses). 
Osteoporosis (OP) and osteoarthritis (OA) occur in a substantial portion of the human 
population over the age of fifty. The National Osteoporosis Foundation (www.nof.org) 
estimates that as many as 44 million Americans are affected by osteoporosis and low 
bone mass. In 1997 the estimated cost for osteoporosis related fractures was $13 billion. 
That figure increased to $17 billion in 2002 and is projected to increase to $210-240 billion 
by 2040. Currently it is expected that one in two women over the age of 50 will suffer an 
osteoporosis-related fracture. 

[003] In spite of its societal impact and prevalence, there is a paucity of 

information on the factors that cause osteoporosis and osteoarthritis to progress more 
rapidly in some individuals than others. Previously considered diseases with little 
opportunity for therapeutic intervention, osteoporosis and osteoarthritis are now 
increasingly viewed as a dynamic process with potential for new pharmacologic and 
surgical treatment modalities. 
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[004] However, the appropriate deployment and selection of treatment 

interventions for OP and OA Is dependent on the development of better methods for the 
assessment of the condition of a patient. Moreover, better methods for the assessment of 
the condition of a patient can also help with drug discovery, drug development, diagnosis, 
prognostication, disease monitoring, therapeutic monitoring and patient management. 

[005] The present invention discloses novel methods and techniques for assessing 

the condition of a bone(s) and/or a joint(s). 

Summary of the Invention 
[006] The invention discloses a method for analyzing at least one of bone mineral 

density, bone structure and surrounding tissue. The method typically comprises: 
(a)obtaining an image of a subject; (b) locating a region of interest on the image; (c) 
obtaining data from the region of interest; and (d) deriving data selected from the group of 
qualitative and quantitative from the image data obtained at step c. 
[007] Further, a kit is provided for aiding in the assessment of the condition of at 

least one of a bone and joint. The kit typically comprises a software program which reads 
at least one of a degeneration pattern, a pattern of normal tissue, a pattern of abnormal 
tissue, and a pattern of diseased tissue. The kit can also include a database of 
measurements for comparison to the at least one of degeneration pattern, pattern of 
normal tissue, pattern of abnormal tissue, and pattern of diseased tissue. Additionally, the 
kit can include a subset of a database of measurements for comparison to the at least one 
of degeneration pattern, pattern of normal tissue, pattern of abnormal tissue, and pattern 
of diseased tissue. 

[008] The Invention also includes at least one of an automated and semi- 

automated method of using an imaging marker. This automatic or semi-automatic method 
comprises: obtaining image data from a subject; obtaining data from the image data 
wherein the data obtained is at least one of quantitative and qualitative data; and 
administering an agent. Either the automated or semi-automated method can be used for: 
drug discovery, diagnosis, disease staging, disease monitoring, disease management, 
prognostication, therapy monitoring, drug efficacy monitoring, and disease prediction. 
[009] In another embodiment, a system for monitoring the efficacy of an agent 

and/or a system for drug discovery is provided. The system includes: administering an 
agent to a subject; obtaining image data; and obtaining data from the image data wherein 
the data obtained is at least one of quantitative and qualitative data. 
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[010] A system is also provided for diagnosing a disease, determining disease 

staging, monitoring disease progression, managing a disease, disease prognostication, 
predicting a disease, monitoring therapy and/or randomizing a subject within a group of 
patients comprising. Any of these systems can include the steps of: (a) obtaining image 
data of a subject; (b) obtaining data from the image data wherein the data obtained is at 
least one of quantitative and qualitative data; and (c) comparing the at least one of 
quantitative and qualitative data in step b to at least one of: a database of at least one of 
quantitative and qualitative data obtained from a group of subjects; at least one of 
quantitative and qualitative data obtained from the subject; and at least one of a 
quantitative and qualitative data obtained from the subject at time Tn. 
[011] In any of these described invention, additional steps can be provided. Such 

additional steps include, for example, enhancing image data. 

[012] Suitable subjects for these steps include for example mammals, humans and 

horses. Suitable anatomical regions of subjects include, for example, dental, spine, hip, 
knee and bone core x-rays. 

[013] A variety of systems can be employed to practice the inventions. Typically at 

least one of the steps of any of the methods is performed on a first computer. Although, it 
is possible to have an arrangement where at least one of the steps of the method is 
performed on a first computer and at least one of the steps of the method is performed on 
a second computer. In this scenario the first computer and the second computer are 
typically connected. Suitable connections include, for example, a peer to peer network, 
direct link, intranet, and internet. 

[014] It is important to note that any or all of the steps of the inventions disclosed 

can be repeated one or more times in series or in parallel with or without the repetition of 
other steps in the various methods. This includes, for example repeating the step of 
locating a region of interest, or obtaining image data. 

[015] Data can also be converted from 2D to 3D to 4D and back; or from 2D to 4D. 

Data conversion can occur at multiple points of processing the information. For example, 
data conversion can occur before or after pattern evaluation and/or analysis. 
[016] Any of the processes described herein are suitable for including the step of 

administering a candidate agent, where that step has not already been performed. 
Suitable agents include, for example, substances administered to a subject, substances 
ingested by a subject, molecules, pharmaceuticals, biopharmaceuticals, 
agropharmaceuticals, genetically manufactured and altered substances, etc. 
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[017] Any data obtained, extracted or generated under any of the methods can be 

compared to a database, a subset of a database, or data previously obtained, extracted or 
generated from the subject. 

[018] The present invention provides methods that allow for the analysis of bone 

mineral density, bone and/or cartilage structure and morphology and/or surrounding tissue 
from images including electronic images and, accordingly, allows for the evaluation of the 
effect(s) of an agent (or agents) on bone and/or cartilage. It is important to note that an 
effect on bone and/or cartilage can occur in agents intended to have an effect, such as a 
therapeutic effect, on bone and/or cartilage as well as agents intended to primarily effect 
other tissues in the body but which have a secondary, or tangential, effect on bone and/or 
cartilage. The images (e.g., x-ray images) can be, for example, dental, hip, spine or other 
radiographs and can be taken from any mammal. The images can be in electronic format. 

[019] The invention includes a method to derive quantitative information on bone 

structure and/or bone mineral density from an image comprising (a) obtaining an image, 
wherein the image optionally includes an external standard for determining bone density 
and/or structure; and (b) analyzing the image obtained in step (a) to derive quantitative 
information on bone structure. The image is taken of a region of interest (ROI). Suitable 
ROI include, for example, a hip radiograph or a dental x-ray obtained on dental x-ray film, 
including the mandible, maxilla or one or more teeth. In certain embodiments, the image 
is obtained digitally, for example using a selenium detector system, a silicon detector 
system or a computed radiography system. In other embodiments, the image can be 
digitized from film, or another suitable source, for analysis. 

[020] A method is included where one or more candidate agents can be tested 

for its effects on bone. Again, the effect can be a primary effect or a secondary effect. For 
example, the candidate agent can be administered to a subject; thereafter an electronic 
image of a portion of a bone of the subject can be obtained; and finally the image 
obtained can be analyzed for Information on bone structure. Information on bone structure 
can relate to a variety of parameters, Including the parameters shown In Table 1 , Table 2 
and Table 3, infra. The images or data can then be compared to a database of images or 
data (e.g., "normal" images or data) and/or compared to one or more images or data taken 
from the same subject prior to administration of the candidate agent or to one or more 
images or data obtained in a reference population. The candidate agent can, for example, 
be molecules, proteins, peptides, naturally occurring substances, chemically synthesized 
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substances, or combinations and cocl<talls thereof. Typically, an agent is one or more 
drugs. Further, the agent can be evaluated for the ability to effect bone diseases such as 
the risk of bone fracture (e.g., osteoporotic fracture). 

[0211 In any of the methods described herein, the analysis can comprise using one 

or more computer programs (or units). Additionally, the analysis can comprise identifying 
one or more regions of interest (ROI) In the image, either prior to, concurrently or after 
analyzing the image, e.g. for information on bone mineral density and/or bone structure. 
The bone density information can be, for example, areas of highest, lowest or median 
density. Bone structural information can be, for example, one or more of the parameters 
shown in Table 1 , Table 2 and Table 3. The various analyses can be performed 
concurrently or in series. Further, when using two or more indices each of the indices can 
be weighted equally or differently, or combinations thereof where more than two indices 
are employed. Additionally, any of these methods can also include analyzing the image 
for bone mineral density information using any of the methods described herein. 

[022] Any of the methods described herein can further comprise applying one or 

more correction factors to the data obtained from the image. For example, correction 
factors can be programmed into a computer unit. The computer unit can be the same one 
that perfonms the analysis of the image or can be a different unit. In certain embodiments, 
the correction factors account for the variation in soft-tissue thicl<ness in individual 
subjects. 

[023] These and other embodiments of the subject invention will readily occur to 

those of skill in the art in light of the disclosure herein. 

Brief Description of the Drawings 

[024] Figs. 1a and b are block diagrams showing the steps for extracting data from 

an image and then deriving quantitative and/or qualitative data from the image. 

[025] Figs. 2a-c are diagrams showing an image taken of a region of anatomical 

interest further illustrating possible locations of regions of interest for analysis. 

[026] Figs. 3a-j illustrate various abnormalities that might occur including, for 

example, cartilage defects, bone marrow edema, subchondral sclerosis, osteophytes and 
cysts. 
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[027] Figs. 4a and b are block diagrams of the method of Fig. 1a showing that the 

steps can be repeated. 

[028] Figs. 5a-e are block diagrams illustrating steps involved in evaluating 

patterns in an image of a region of interest. 

[029] Fig. 6a-e are block diagrams Illustrating steps involved in deriving 

quantitative and qualitative data from an Image in conjunction with administering 
candidate molecules or drugs for evaluation, 

[030] Figs. 7a-d are block diagrams illustrating steps involved in comparing 

derived quantitative and qualitative information to a database or to information obtained at 
a previous time. 

[031] Figs. 8a-d are block diagrams illustrating steps involved in comparing 

converting an image to a pattern of normal and diseased tissue 

[032] Fig, 9 is a diagram showing the use one or more devices in the process of 

developing a degeneration pattern and using a database for degeneration patterns. 

Detailed Description 

[033] The following description is presented to enable any person skilled in the art 

to make and use the invention. Various modifications to the embodiments described will 
be readily apparent to those skilled in the art, and the generic principles defined herein 
can be applied to other embodiments and applications without departing from the spirit 
and scope of the present invention as defined by the appended claims. Thus, the present 
invention is not intended to be limited to the embodiments shown, but is to be accorded 
the widest scope consistent with the principles and features disclosed herein. To the 
extent necessary to achieve a complete understanding of the invention disclosed, the 
specification and drawings of all issued patents, patent publications, and patent 
applications cited in this application are incorporated herein by reference, 

[034] The practice of the present invention employs, unless otherwise indicated, 

currently conventional methods of imaging and image processing within the skill of the art. 
Such techniques are explained fully in the literature. See, e.g., WO 02/22014, X-Ray 
Structure Determination: A Practical Guide, 2nd Edition, editors Stout and Jensen, 1989, 
John Wiley & Sons, publisher; Body CT: A Practical Approach, editor Slone, 1999, 
McGraw-Hill publisher; The Essential Physics of Medical Imaging, editors Bushberg, 
Seibert, Leidholdt Jr & Boone, 2002, Lippincott, Williams & Wilkins; X-ray Diagnosis: A 
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Physician's Approach, editor Lam, 1998 Springer-Verlag, publisher; Dental Radiology: 
Understanding the X-Ray Image, editor Laetitia Brocklebank 1997, Oxford University 
Press publisher; and Digital Image Processing, editor Kenneth R. Castleman, 1996 
Prentice Hall, publisher; The Image Processing Handbook, editor John C. Russ, 3rd 
Edition, 1998, CRC Press; Active Contours: The Application of Techniques from Graphics, 
Vision, Control Theory and Statistics to Visual Tracking of Shapes in Motion, Editors 
Andrew Blake, Michael Isard, 1 999 Springer Verlag. As will be appreciated by those of skill 
in the art, as the field of imaging continues to advance methods of imaging currently 
employed can evolve over time. Thus, any imaging method or technique that is currently 
employed is appropriate for application of the teachings of this invention as'well as 
techniques that can be developed in the future. A further detailed description of imaging 
methods is not provided in order to avoid obscuring the invention. 

[035] As shown in Fig. 1a, the first step is to locate a part of the body of a subject, 

for example in a human body, for study 98. The part of the body located for study is the 
region of anatomical interest (RAI). In locating a part of the body for study, a determination 
Is made to, for example, take an image or a series of Images of the body at a particular 
location, e.g. hip, dental, spine, etc. Images include, for example, conventional x-ray 
images, x-ray tomosynthesis, ultrasound (including A-scan, B-scan and C-scan) computed 
tomography (CT scan), magnetic resonance imaging (MRl), optical coherence 
tomography, single photon emission tomography (SPECT), and positron emission 
tomography, or such other imaging tools that a person of skill in the art would find useful in 
practicing the invention. Once the image is taken, a region of interest (ROI) can be 
located within the image 100. Image data is extracted from the image 102. Finally, 
quantitative and/or qualitative data is extracted from the image data 120, The quantitative 
and/or qualitative data extracted from the image includes, for example, the parameters 
and measurements shown in Table 1, Table 2 or Table 3. 

[036] Each step of locating a part of the body for study 98, optionally locating a 

region of interest 100, obtaining image data 102, and deriving data 120, can be repeated 
one or more times 99,101, 103, 121, respectively, as desired. 

[037] As shown in Fig. 1b image data can be optionally enhanced 104 by applying 

image processing techniques, such as noise filtering or diffusion filtering, to facilitate 
further analysis. Similar to the process shown in Fig. 1a, locating a part of the body for 
study 98, optionally locating a region of interest 100, obtaining image data 102, enhancing 
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image data 104, and deriving data 120, can be repeated one or more times 99,101, 103, 
105, 121, respectively, as desired. 



Table 1 

Representative Parameters Measured with 
Quantitative and Qualitative Image Analysis Methods 



Parameter 


Measurements 


DUl ic {Jciiciiiiddo 


* 


oiainieos sieei ecjuivoicni iniuKness 






^Average gray vaiue oi ine region or inieresi expressea as inicKness 






OT siainisss sieei wiin squivaieni intensiiyj 




• 


1 raDecuiar conirasi 




• 


(Trabecular equivalent thickness / Marrow equivalent thickness) 




• 


Fractal dimension 




• 


Fourier spectral analysis 






(Mean transform coefficient absolute value and mean spatial first 






moment) 




• 


Predominant orientation of spatial energy spectrum 




• 


Trabecular area 






(Pixel count of extracted trabeculae) 




• 


Trabecular area / Total area 




• 


Trabecular perimeter 






(Count of trabecular pixels with marrow pixels in their neighborhood, 






proximity or vicinity) 




• 


Trabecular distance transform 






(For each trabecular pixel, calculation of distance to closest marrow 






pixel) 




• 


Marrow distance transform 






(For each marrow pixel, calculation of distance to closest trabecular 






pixel) 




• 


Trabecular distance transform regional maximal values (mean, min., 






max, std. Dev). 
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directions) 




• 


Trabecular Bone Pattern Factor 






(TBPf = (P1 - P2) / (A1 - A2 ) where P1 and A1 are the perimeter 






length and trabecular bone area before dilation and P2 and A2 






corresponding values after a single pixel dilation, measure of 






connectivity) 






Connected skeleton count or Trees (T) 






Node count (N) 






Segment count (S) 






Node-to-node segment count (NN) 






Node-to-free-end segment count (NF) 






Node-to-node segment length (NNL) 






Node-to-free-end segment length (NFL) 






Free-end-to-free-end segment length (FFL) 






Node-to-node total struts length (NN.TSL) 






Free-end-to-free-ends total struts length( FF.TSL) 






Total struts length (TSL) 






FF.TSL/ TSL 






NN.TSL/ TSL 
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Parameter 


Measurements 




• Loop count (Lo) 

• Loop area 

• Mean distance transform vaiues for each connected skeleton 

• Mean distance transform values for each segment (Tb.Th ) 

• Mean distance transform values for each node-to-node segment 
(Tb.Th. NN) 

• Mean distance transform values for each node-to-free-end segment 
(Tb.Th. NF) 

• Orientation (angle) of each segment 

• Angle between segments 

• Length-thickness ratios (NNL/Tb.Th.NN ) and (NFL/ Tb.Th.NF) 

• Interconnectivity index (ICI) ICI = (N * NN)/ ( T * (NF + 1 ) ) 


Cartilage and 
cartilage 
defect/diseased 
cartilage parameters 


• Totai cartilage volume 

• Partial/Focal cartilage volume 

• Cartilage thickness distribution (thickness map) 

• Mean cartilage thickness for total region or focal region 

• Median cartilage thickness for totai region or focal region 

• Maximum cartilage thickness for total region or focal region 

• Minimum cartilage thickness for total region or focal region 

• 3D cartilage surface Information for total region or focal region 

• Cartilage curvature analysis for total region or focal region 

• Volume of cartilage defect/diseased cartilage 

• Depth of cartilage defect/diseased cartilage 

• Area of cartilage defect/diseased cartilage 

• 2D or 3D location of cartilage defect/diseased cartilage in articular 
surface 

• 2D or 3D location of cartilage defect/diseased cartilage In 
relationship to weight-bearing area 

• Ratio: diameter of cartilage defect or diseased cartilage / 
thickness of surrounding normal cartilage 

• Ratio: depth of cartilage defect or diseased cartilage / thickness 
of surrounding normal cartilage 

• Ratio: volume of cartilage defect or diseased cartilage / 
thickness of surrounding normal cartilage 

• Ratio: surface area of cartilage defect or diseased cartilage / 
total joint or articular surface area 

• Ratio: volume of cartilage defect or diseased cartilage / total 
cartilage volum^ 


Other articuiar 
parameters 


• Presence or absence of bone marrow edema 

• Volume of bone marrow edema 

• Volume of bone marrow edema normalized by width, area, size, 
volume of femoral condyle(s)/tlbial plateau/patella - other bones In other 
joints 

• Presence or absence of osteophytes 

• Presence or absence of subchondral cysts 

• Presence or absence of subchondral sclerosis 

• Volume of osteophytes 

• Volume of subchondral cysts 

• Volume of subchondral sclerosis 

• Area of bone marrow edema 

• Area of osteophytes 

• Area of subchondral cysts 

• Area of subchondral sclerosis 

• Depth of bone marrow edema 

• Depth of osteophytes 

• Depth of subchondral cysts 
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Parameter 


Measurements 




• Depth of subchondral sclerosis 

• Volume, area, depth of osteophytes, subchondral cysts, subchondral 
sclerosis normalized by width, area, size, volume of femoral 
condyle(s)/tibial plateau/patella - other bones in other joints 

• Presence or absence of meniscal tear 

• Presence or absence of cruciate ligament tear 

• Presence or absence of collateral ligament tear 

• Volume of menisci 

• Ratio of volume of normal to torn/damaged or degenerated meniscal 
tissue 

• Ratio of surface area of normal to torn/damaged or degenerated 
meniscal tissue 

• Ratio of surface area of normal to torn/damaged or degenerated 
meniscal tissue to total joint or cartilage surface area 

• Ratio of surface area of torn/damaged or degenerated meniscal 
tissue to total joint or cartilage surface area 

• Size ratio of opposing articular surfaces 

• Meniscal subluxation/dislocation in mm 

• Index combining different articular parameters which can also 
include 

o Presence or absence of cruciate or collateral ligament 
tear 

o Body mass Index, weight, height 

• 3D surface contour information of subchondral bone 

• Actual or predicted knee flexion angle during gait cycle 
(latter based on gait patterns from subjects with matching 
demographic data retrieved from motion profile database) 

• Predicted knee rotation during gait cycle 

• Predicted knee displacement during gait cycle 

• Predicted load bearing line on cartilage surface during gait cycle and 
measurement of distance between load bearing line and cartilage 
defect/diseased cartilage 

• Predicted load bearing area on cartilage surface during gait cycle 
and measurement of distance between load bearing area and 
cartilage defect/diseased cartilage 

• Predicted load bearing line on cartilage surface during standing or 
different degrees of knee flexion and extension and measurement of 
distance between load bearing line and cartilage defect/diseased 
cartilage 

• Predicted load bearing area on cartilage surface during standing or 
different degrees of knee flexion and extension and measurement of 
distance between load bearing area and cartilage defect/diseased 
cartilage 

• Ratio of load bearing area to area of cartilage defect/diseased 

cartilage 

• Percentage of load bearing area affected by cartilage disease 

• Location of cartilage defect within load bearing area 

• Load applied to cartilage defect, area of diseased cartilage 

• Load applied to cartilage adjacent to cartilage defect, area of 
diseased cartilage 



[038] As will be appreciated by tliose of skill in tlie art, tine parameters and 

nneasurements shown in Table 1 are provided for illustration purposes. Other parameters 
and measurements, ratios, derived values or indices can be used to extract quantitative 
and/or qualitative information about the ROI without departing from the scope of the 
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invention. Additionally, where multiple ROI or multiple derivatives of data are used, the 
parameter measured can be the same parameter or a different parameter without 
departing from the scope of the invention. Additionally, data from different ROI's can be 
combined or compared as desired. 

[039] Additional measurements can be performed that are selected based on the 

anatomical structure to be studied as described below. 

[040] Once the data is extracted from the Image it can be manipulated to assess 

the severity of the disease and to determine disease staging (e.g., mild, moderate, severe 
or a numerical value or Index). The information can also be used to monitor progression of 
the disease and/or the efficacy of any interventional steps that have been taken. Finally, 
the information can be used to predict the progression of the disease or to randomize 
patient groups in clinical trials. 

[041] Fig. 2a illustrates an image 200 taken of an RAI, shown as 202. As shown in 

Fig. 2a, a single region of interest (ROI) 210 has been identified within the image. The 
ROI 210 can take up the entire image 200, or nearly the entire image. As shown in 
Fig. 2b more than one ROI can be identified in an image. In this example, a first ROI 220 
is depicted in one region of the image 200, and a second ROI 222 is depicted within the 
image. In this instance, neither of these ROI overlap or abut As will be appreciated by a 
person of skill in the art, the number of ROI identified in an image 200 is not limited to the 
two depicted. Turning now to Fig. 2c another embodiment showing two ROI for illustration 
purposes is shown. In this instance, the first ROI 230 and the second ROI 232, are 
partially overlapping. As will be appreciated by those of skill in the art, where multiple ROI 
are used any or all of the ROI can be organized such that it does not overlap, it abuts 
without overlapping, it overlaps partially, it overlaps completely (for example where a first 
ROI is located completely within a second identified ROI), and combinations thereof. 
Further the number of ROI per image 200 can range from one (ROIi) to n (ROIn) where n 
is the number of ROI to be analyzed. 

[042] Bone density, microarchitecture, macro-anatomic and/or biomechanical (e.g. 

derived using finite element modeling) analyses can be applied within a region of 
predefined size and shape and position. This region of interest can also be referred to as 
a "window." Processing can be applied repeatedly within the window at different positions 
of the image. For example, a field of sampling points can be generated and the analysis 
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performed at these points. The results of the analyses for each paranneter can be stored 
in a matrix space, e.g., where its position corresponds to the position of the sampling point 
where the analysis occurred, thereby forming a map of the spatial distribution of the 
parameter (a parameter map). The sampling field can have regular intervals or irregular 
intervals with varying density across the image. The window can have variable size and 
shape, for example to account for different patient size or anatomy. 

[043] The amount of overlap between the windows can be determined, for 

example, using the interval or density of the sampling points (and resolution of the 
parameter maps). Thus, the density of sampling points is set higher in regions where 
higher resolution is desired and set lower where moderate resolution is sufficient, in order 
to improve processing efficiency. The size and shape of the window would determine the 
local specificity of the parameter. Window size is preferably set such that it encloses most 
of the structure being measured. Oversized windows are generally avoided to help 
ensure that local specificity is not lost. 

[044] The shape of the window can be varied to have the same orientation and/or 

geometry of the local structure being measured to minimize the amount of structure 
clipping and to maximize local specificity. Thus, both 2D and/or 3D windows can be used, 
as well as combinations thereof, depending on the nature of the image and data to be 
acquired. 

[045] In another embodiment, bone density, microarchitecture, macro-anatomic 

and/or biomechanical (e.g. derived using finite element modeling) analyses can be applied 
within a region of predefined size and shape and position. The region is generally 
selected to include most, or all, of the anatomic region under investigation and, preferably, 
the parameters can be assessed on a pixel-by-pixel basis {e.g., in the case of 2D or 3D 
images) or a voxel-by-voxel basis in the case of cross-sectional or volumetric images 
(e.g., 3D images obtained using MR and/or CT). Alternatively, the analysis can be applied 
to clusters of pixels or voxels wherein the size of the clusters is typically selected to 
represent a compromise between spatial resolution and processing speed. Each type of 
analysis can yield a parameter map. 

[046] Parameter maps can be based on measurement of one or more parameters 

in the image or window; however, parameter maps can also be derived using statistical 
methods. In one embodiment, such statistical comparisons can include comparison of 
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data to a reference population, e.g. using a z-score or a T-score. Thus, parameter maps 
can include a display of z-scores or T-scores. 

[047] Additional measurements relating to the site to be measured can also be 

taken. For example, measurements can be directed to dental, spine, hip, knee or bone 
cores. Examples of suitable site specific measurements are shown in Table 2. 



Table 2 

Site specific measurement of bone parameters 



Parameters sDocific to 
hiD x-ravs 


• All microarchitecture parameters on structures parallel to 
stress lines 

• All microarchitecture parameters on structures perpendicular 
to stress lines 

Geometry 

• Shaft angle 

• Neck angle 

• Diameter of femur neck 

• Hip axis length 

• Largest cross-section of femur head 

• Average thickness of cortical within ROI 

• Standard deviation of cortical thickness within ROI 

• IVlaximum thickness of cortical within ROI 

• Minimum thickness of cortical within ROI 
Hip joint space width 


Parameters specific to 
SDine x-ravs 


• All microarchitecture parameters on vertical structures 

• All microarchitecture parameters on horizontal structures 
Geometry 

• Superior endplate cortical thickness (anterior, center, 

pUoLc;! lUI ) 

• Inferior endplate cortical thickness (anterior, center, 

posterior) 

• Anterior vertebral wall cortical thickness (superior, center, 
inferior) 

• Posterior vertebral wall cortical thickness (superior, center, 
inferior) 

• Superior aspect of pedicle cortical thickness 

• inferior aspect of pedicle cortical thickness 

• Vertebral height (anterior, center, posterior) 

• Vertebral diameter (superior, center, Inferior), 

• Pedicle thickness (supero-inferior direction). 

• Maximum vertebral height 

• Minimum vertebral height 

• Average vertebral height 

• Anterior vertebral height 

• Medial vertebral height 

• Posterior vertebral height 

• Maximum inter-vertebra! height 

• Minimum Inter-vertebral height 

• Average inter-vertebral height 


Parameters SDeclfic to 
knee x-ravs 


• Average medial joint space width 

• Minimum medial joint space width 

• Maximum medial joint space width 

• Average lateral joint space width 

• Minimum lateral joint space width 

• Maximum lateral Joint space width 
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[048] As will be appreciated by those of skill in the art, measurement and Image 

processing techniques are adaptable to be applicable to both microarchitecture and 
macro-anatomical structures. Examples of these measurements are shown in Table 3. 
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Table 3 

Measurements applicable on Microarchitecture and Macro-anatomical Structures 



Averaae densltv 
measurement 


• Calibrated density of ROl 


Measurements on 
micro-anatomical 
structures of dental, 
spine, hiD. knee or bone 
cores x-ravs 


The following parameters are derived from the extracted structures: 

• Calibrated density of extracted structures 

• Calibrated density of background 

• Average Intensity of extracted structures 

• Average intensity of background (area other than extracted 
structures) 

• Structural" contrast (average intensity of extracted structures 
/ average intensity of background ) 

• Calibrated structural contrast (calibrated density extracted 
structures / calibrated density of background) 

• Total area of extracted structures 

• Total area of ROl 

• Area of extracted structures normalized by total area of ROl 

• Boundary lengths (perimeter) of extracted normalized by 

total flr^a nf ROl 

• Number of structures normalized by area of ROl 

• Trabecular bone pattern factor; measures concavity and 
convexity of structures 

• Star volume of extracted structures 

• Star volume of background 

• Number of loops normalized by area of ROl 


Measurements on 
Distance transform 
of extracted 
structures 


The following statistics are measured from the distance transform 
regional maximum values: 

• Average regional maximum thickness 

• Standard deviation of regional maximum thickness 

• Largest value of regional maximum thickness 

• Median of regional maximum thickness 
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Measurements on 
skeleton of 
extracted structures 


• Average length of networks (units of connected segments) 

• Maximum length of networks 

• Average thickness of structure units (average distance 
transform values along skeleton) 

• Maximum thickness of structure units (maximum distance 
transform values along skeleton) 

• Number of nodes normalized by ROI area 

• Number of segments normalized by ROI area 

• Number of free-end segments normalized by ROI area 

• Number of inner (node-to-node) segments normalized ROI 
area 

• Average segment lengths 

• Average free-end segment lengths 

• Average Inner segment lengths 

• Average orientation angle of segments 

• Average orientation angle of inner segments 

• Segment tortuosity; a measure of straightness 

• Segment solidity; another measure of straightness 

• Average thickness of segments (average distance transform 
values along skeleton segments) 

• Average thickness of free-end segments 

• Average thickness of inner segments 

• Ratio of inner segment lengths to inner segment thickness 

• Ratio of free-end segment lengths to free-end segment 
thickness 

• Interconnectivity index; a function of number of inner 
segments, free-end segments and number of networks. 


Directional siceleton 

segment 

measurements 


All measurement of skeleton segments can be constrained by 
one or more desired orientation by measuring only skeleton 
segments within ranges of angle. 


Waterslied 
segmentation 


Watershed segmentation is applied to gray level images. 
Statistics of watershed segments are; 

• Total area of segments 

• Number of segments normalized by total area of segments 

• Average area of segments 

• Standard deviation of segment area 

• Smallest segment area 

• Largest segment area 



[049] Calibrated density typically refers to the measurement of intensity values of 

features in Images converted to its actual material density or expressed as the density of a 
reference material whose density is known. The reference material can be metal, polymer, 
plastics, bone, cartilage, etc., and can be part of the object being imaged or a calibration 
phantom placed in the imaging field of view during image acquisition. 

[050] Extracted structures typically refer to simplified or amplified representations 

of features derived from images. An example would be binary images of trabecular 
patterns generated by background subtraction and thresholding. Another example would 
be binary images of cortical bone generated by applying an edge filter and thresholding. 
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The binary images can be superimposed on gray level images to generate gray level 
patterns of structure of interest. 

[051] Distance transform typically refers to an operation applied on binary images 

where maps representing distances of each 0 pixel to the nearest 1 pixel are generated. 
Distances can be calculated by the Euclidian magnitude, city-block distance, La Place 
distance or chessboard distance. 

[052] Distance transform of extracted structures typically refer to distance 

transform operation applied to the binary images of extracted structures, such as those 
discussed above with respect to calibrated density. 

[053] Skeleton of extracted structures typically refer to a binary image of 1 pixel 

wide patterns, representing the centeriine of extracted structures. It Is generated by 
applying a skeletonization or medial transform operation, by mathematical morphology or 
other methods, on an image of extracted structures. 

[054] Skeleton segments typically are derived from skeleton of extracted 

structures by performing pixel neighborhood analysis on each skeleton pixel. This analysis 
classifies each skeleton pixel as a node pixel or a skeleton segment pixel. A node pixel 
has more than 2 pixels in its 8-neighborhood. A skeleton segment is a chain of skeleton 
segment pixels continuously 8-connected. Two skeleton segments are separated by at 
least one node pixel. 

[055] Watershed segmentation as it is commonly known to a person of skill in the 

art, typically is applied to gray level images to characterize gray level continuity of a 
structure of interest. The statistics of dimensions of segments generated by the process 
are, for example, those listed in Table 3 above. As will be appreciated by those of skill in 
the art, however, other processes can be used without departing from the scope of the 
invention. 

[056] Turning now to Fig. 3a, a cross-section of a cartilage defect is shown 300. 

The cross-hatched zone 302 corresponds to an area where there is cartilage loss. FiQ. 3b 
is a top view of the cartilage defect shown in Fig. 3a. 

[057] Fig. 3c illustrates the depth of a cartilage defect 310 in a first cross-section 

dimension with a dashed line illustrating a projected location of the original cartilage 
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surface 312. By comparing these two values a ratio of cartilage defect deptli to cartilage 
defect width can be calculated. 

[058] Fig.Sd illustrated the depth of the cartilage 320 along with the width of the 

cartilage defect 322. These two values can be compared to determine a ratio of cartilage 
depth to cartilage defect width. 

[059] Fig. 3e shows the depth of the cartilage defect 310 along with the depth of 

the cartilage 320. A dashed line is provided illustrating a projected location for the original 
cartilage surface 3f 2. Similar to the measurements made above, ratios between the 
various measurements can be calculated. 

[060] Turning now to Fig. 3f, an area of bone marrow edema is shown on the 

femur 330 and the tibia 332. The shaded area of edema can be measured on a 
T2-weighted MRI scan. Alternatively, the area can be measured on one or more slices. 
These measurements can then be extended along the entire joint using multiple slices or 
a 3D acquisition. From these measurements volume can be determined or derived, 

[061] Fig. 3g shows an area of subchondral sclerosis In the acetabulum 340 and 

the femur 342. The sclerosis can be measured on, for example, a T1 or T2-weighted MRI 
scan or on a CT scan. The area can be measured on one or more slices. Thereafter the 
measurement can be extended along the entire joint using multiple slices or a 3D 
acquisition. From these values a volume can be derived of the subchondral sclerosis. For 
purposes of illustration, a single sclerosis has been shown on each surface. However, a 
person of skill in the art will appreciate that more than one sclerosis can occur on a single 
joint surface, 

[062] Fig,3h shows osteophytes on the femur 350 and the tibia 352. The 

osteophytes are shown as cross-hatched areas. Similar to the sclerosis shown in Fig. 3g, 
the osteophytes can be measured on, for example, a T1 or T2-weighted MRI scan or on a 
CT scan. The area can be measured on one or more slices. Thereafter the measurement 
can be extended along the entire joint using multiple slices or a 3D acquisition. From 
these values a volume can be derived of the osteophytes. Additionally, a single 
osteophyte 354 or osteophyte groups 356 can be included in any measurement. Persons 
of skill in the art will appreciate that groups can be taken from a single joint surface or from 
opposing joint surfaces, as shown, without departing from the scope of the invention. 
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[063] Turning now to Fig. 3i an area of subchondral cysts 360, 362, 364 is sinown. 

Similar to the sclerosis shown In Fig. 3g, the cysts can be measured on, for example, a T1 
or T2-welghted MRI scan or on a CT scan. The area can be measured on one or more 
slices. Thereafter the measurement can be extended along the entire joint using multiple 
slices or a 3D acquisition. From these values a volume can be derived of the cysts. 
Additionally, single cysts 366 or groups of cysts 366' can be included in any 
measurement. Persons of skill in the art will appreciate that groups can be taken from a 
single joint surface, as shown, or from opposing joint surfaces without departing from the 
scope of the invention. 

[064] Fig. 3j illustrates an area of torn meniscal tissue (cross-hatched) 372, 374 

as seen from the top 370 and in cross-section 371. Again, similar to the sclerosis shown 
in Fig. 3g, the torn meniscal tissue can be measured on, for example, a T1 or T2-weighted 
MRI scan or on a CT scan. The area can be measured on one or more slices. Thereafter 
the measurement can be extended along the entire joint using multiple slices or a 3D 
acquisition. From these values a volume can be derived of the tear. Ratios such as 
surface or volume of torn to normal meniscal tissue can be derived as well as ratios of 
surface of torn meniscus to surface of opposing articulating surface. 

[065] As shown in Fig. 4a, the process of optionally locating a ROI 100, extracting 

image data from the RO1 102, and deriving quantitative and/or qualitative image data from 
the extracted image data 120, can be repeated 122. Alternatively, or in addition, the 
process of locating a ROI 100, can be repeated 124, A person of skill in the art will 
appreciate that these steps can be repeated one or more times in any appropriate 
sequence, as desired, to obtain a sufficient amount of quantitative and/or qualitative data 
on the ROI or to separately extract or evaluate parameters. Further, the ROI used can be 
the same ROI as used in the first process or a newly identified ROI In the image. 
Additionally, as with Fig. 1a the steps of locating a region of interest 100, obtaining image 
data 102, and deriving quantitative and/or qualitative image data can be repeated one or 
more times, as desired, 101, 103, 121, respectively. Although not depicted here, as 
discussed above with respect to Fig. 1a, the additional step of locating a part of the body 
for study 98 can be performed prior to locating a region of Interest 100 without departing 
from the invention. Additionally that step can be repeated 99. 

[066] Fig. 4b illustrates the process shown in Fig. 4a with the additional step 

enhancing image data 104. Additionally, the step of enhancing image data 104 can be 
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repeated one or more times 105, as desired. The process of enliancing image data 104 
can be repeated 126 one or more times as desired. 

[067] Turning now to Fig. 5a, a process is sliown wliereby a region of interest is 

optionally located 100. Althougli not depicted liere, as discussed above witln respect to 
Fig. 1a, the step of locating a part of the body for study 98 can be performed prior to 
locating a region of interest ^00 without departing from the invention. Additionally that step 
can be repeated 99. Once the region of interest is located 100, and image data is 
extracted from the ROI 102, the extracted image data can then be converted to a 2D 
pattern 130, a 3D pattern 132 or a 4D pattern 133, for example including velocity or time, 
to facilitate data analyses. Following conversion to 2D 130, 3D 132 or 4D pattern 133 the 
images are evaluated for patterns 140, Additionally images can be converted from 2D to 
3D 131, or from 3D to 4D 131', if desired. Although not illustrated to avoid obscuring the 
figure, persons of skill in the art will appreciate that similar conversions can occur between 
2D and 4D in this process or any process illustrated in this invention. 

[068] As will be appreciated by those of skill in the art, the conversion step is 

optional and the process can proceed directly from extracting image data from the ROI 
102 to evaluating the data pattern 140 directly 134. Evaluating the data for patterns, 
includes, for example, performing the measurements described in Table 1 ,Table 2 or 
Table 3, above. 

[069] Additionally, the steps of locating the region of interest 100, obtaining image 

data 102, and evaluating patterns 141 can be performed once or a plurality of times, 101, 
103, 141, respectively at any stage of the process. As will be appreciated by those of skill 
in the art, the steps can be repeated. For example, following an evaluation of patterns 
140, additional image data can be obtained 135, or another region of interest can be 
located 137. These steps can be repeated as often as desired, in any combination 
desirable to achieve the data analysis desired. 

[070] Fig. 5b illustrates an alternative process to that shown In Fig. 5a which 

includes the step of enhancing image data 104 prior to converting an image or image data 
to a 2D 130, 3D 132, or 4D 133 pattern. The process of enhancing Image data 104, can 
be repeated 105 if desired. Fig. 5c Illustrates an alternative embodiment to the process 
shown in Fig. 5b. In this process, the step of enhancing image data 104 occurs after 
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converting an image or image data to a 2D 130, 3D 132, or 4D 133 pattern. Again, the 
process of enhancing image data 104, can be repeated 105 if desired. 

[071] Fig. 5d illustrates an alternative process to that shown in Fig. 5a. After 

locating a part of the body for study 98 and Imaging, the image is then converted to a 2D 
pattern 130, 3D pattern 132 or 4D pattern 133, The region of interest 100 is optionally 
located within the Image after conversion to a 2D, 3D or 4D image and data is then 
extracted 102, Patterns are then evaluated in the extracted Image data 140. As with the 
process of Fig. 5a, the conversion step is optional. Further, if desired, images can be 
converted between 2D, 3D 131 and 4D f3f 'if desired. 

[072] Similar to Fig. 5a, some or all the processes can be repeated one or more 

times as desired. For example, locating a part of the body for study 98, locating a region of 
interest 100, obtaining image data 102, and evaluating patterns 140, can be repeated one 
or more times if desired, 99, 101, 103, 141, respectively. Again steps can be repeated. 
For example, following an evaluation of patterns 140, additional Image data can be 
obtained 135, or another region of interest can be located 137 and/or another portion of 
the body can be located for study 139, These steps can be repeated as often as desired, 
in any combination desirable to achieve the data analysis desired. 

[073] Fig. 5e illustrates an alternative process to that shown in Fig. 5d. In this 

process image data can be enhanced 104. The step of enhancing image data can occur 
prior to conversion 143, prior to locating a region of interest 145, prior to obtaining Image 
data 102, or prior to evaluating patterns 149, 

[074] Similar to Fig. 5a, some or all the processes can be repeated one or more 

times as desired, including the process of enhancing image data 104, which is shown as 
105. 

[075] The method also comprises obtaining an image of a bone or a joint, 

optionally converting the image to a two-dimensional or three-dimensional or four- 
dimensional pattern, and evaluating the amount or the degree of normal, diseased or 
abnormal tissue or the degree of degeneration in a region or a volume of interest using 
one or more of the parameters specified in Table 1 , Table 2 and/or Table 3. By performing 
this method at an initial time Ti, information can be derived that is useful for diagnosing 
one or more conditions or for staging, or determining, the severity of a condition. This 
information can also be useful for determining the prognosis of a patient, for example with 
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osteoporosis or arthritis. By performing this method at an Initial time Ti, and a later time 
T2, the change, for example in a region or volume of interest, can be determined which 
then facilitates the evaluation of appropriate steps to take for treatment. IVIoreover, if the 
subject is already receiving therapy or if therapy is Initiated after time Ti, it is possible to 
monitor the efficacy of treatment. By performing the method at subsequent times, T2-Tn. 
additional data ca be acquired that facilitate predicting the progression of the disease as 
well as the efficacy of any interventional steps that have been taken. As will be 
appreciated by those of skill in the art, subsequent measurements can be taken at regular 
time Intervals or irregular time intervals, or combinations thereof. For example, it can be 
desirable to perform the analysis at Ti with an initial follow-up, T2, measurement taken one 
month later. The pattern of one month follow-up measurements could be performed for a 
year (12-one-month intervals) with subsequent follow-ups performed at 6 month Intervals 
and then 12 month intervals. Alternatively, as an example, three initial measurements 
could be at one month, followed by a single six month follow up which is then followed 
again by one or more one month follow-ups prior to commencing 12 month follow ups. 
The combinations of regular and irregular intervals are endless, and are not discussed 
further to avoid obscuring the invention. 

[076] IVIoreover, one or more of the parameters listed in Tables 1 , 2 and 3 can be 

measured. The measurements can be analyzed separately or the data can be combined, 
for example using statistical methods such as linear regression modeling or correlation. 
Actual and predicted measurements can be compared and correlated. 

[077] The method for assessing the condition of a bone or joint in a subject can 

be fully automated such that the measurements of one or more of the parameters 
specified In Table 1, Table 2 or Table 3 are done automatically without intervention. The 
automatic assessment then can include the steps of diagnosis, staging, prognostication or 
monitoring the disease or diseases, or to monitor therapy. As will be appreciated by those 
of skill in the art, the fully automated measurement is, for example, possible with image 
processing techniques such as segmentation and registration. This process can include, 
for example, seed growing, thresholding, atlas and model based segmentation methods, 
live wire approaches, active and/or deformable contour approaches, contour tracking, 
texture based segmentation methods, rigid and non-rigid surface or volume registration, 
for example based on mutual information or other similarity measures. One skilled in the 
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art will readily recognize other techniques and methods for fully automated assessment of 
the parameters and measurements specified in Table 1 , Table 2 and Table 3. 

[078] Alternatively, the method of assessing the condition of a bone or joint in a 

subject can be semi-automated such that the measurements of one or more of the 
parameters, such as those specified in Table 1, are performed semi-automatically, i.e., 
with intervention. The semi-automatic assessment then allows for human interaction and, 
for example, quality control, and utilizing the measurement of said parameter(s) to 
diagnose, stage, prognosticate or monitor a disease or to monitor a therapy. The semi- 
automated measurement is, for example, possible with image processing techniques such 
as segmentation and registration. This can include seed growing, thresholding, atlas and 
model based segmentation methods, live wire approaches, active and/or deformable 
contour approaches, contour tracking, texture based segmentation methods, rigid and 
non-rigid surface or volume registration, for example base on mutual information or other 
similarity measures. One skilled in the art will readily recognize other techniques and 
methods for semi-automated assessment of the parameters specified in Table 1 , Table 2 
or Table 3. 

[079] Turning now to Fig. 6a, a process is shown whereby the user locates a ROI 

100, extracts image data from the ROI 102, and then derives quantitative and/or 
qualitative image data from the extracted image data 120, as shown above with respect to 
Fie. 1. Following the step of deriving quantitative and/or qualitative image data, a 
candidate agent is administered to the patient 150. The candidate agent can be any agent 
the effects of which are to be studied. Agents can include any substance administered or 
ingested by a subject, for example, molecules, pharmaceuticals, biopharmaceuticals, 
agropharmaceuticais, or combinations thereof, including cocktails, that are thought to 
affect the quantitative and/or qualitative parameters that can be measured in a region of 
interest. These agents are not limited to those intended to treat disease that affects the 
musculoskeletal system but this Invention Is Intended to embrace any and all agents 
regardless of the intended treatment site. Thus, appropriate agents are any agents 
whereby an effect can be detected via imaging. The steps of locating a region of interest 
100, obtaining image data 102, obtaining quantitative and/or qualitative data from image 
data 120, and administering a candidate agent 150, can be repeated one or more times as 
desired, 101, 103, 121, 151, respectively. 
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[080] Fig. 6b shows the additional step of enhancing image data 104, which can 

also be optionally repeated 105 as often as desired. 

[081] As shown in Fig, 6c these steps can be repeated one or more times 152 to 

determine the effect of the candidate agent. As will be appreciated by those of skill in the 
art, the step of repeating can occur at the stage of locating a region of interest 152 as 
shown in Fig. 6b or it can occur at the stage obtaining image data 153 or obtaining 
quantitative and/or qualitative data from image data 154 as shown in Fig. 6d. 

[082] Fig. 6e shows the additional step of enhancing image data 104, which can 

optionally be repeated 105, as desired, 

[083] As previously described, some or all the processes shown In Figs. 6a-e can 

be repeated one or more times as desired. For example, locating a region of interest 100, 
obtaining Image data 102, enhancing image data 104, obtaining quantitative and/or 
qualitative data 120, evaluating patterns 140, and administering candidate agent ^50 can 
be repeated one or more times if desired, 101, 103, 105, 121, 141, 151 respectively. 

[084] In the scenario described in relation to Figs. 6, an image is taken prior to 

administering the candidate agent. However, as will be appreciated by those of skill in the 
art, it is not always possible to have an image prior to administering the candidate agent. 
In those situations, progress is determined over time by evaluating the change in 
parameters from extracted image to extracted image. 

[085] Turning now to Fig. 7a, the process is shown whereby the candidate agent 

is administered first 150. Thereafter a region of interest is located in an image taken 100 
and image data is extracted 102, Once the image data is extracted, quantitative and/or 
qualitative data is extracted from the image data 120. In this scenario, because the 
candidate agent is administered first, the derived quantitative and/or qualitative data 
derived is compared to a database 160 or a subset of the database which includes data 
for subjects having similar tracked parameters. As shown in Fig. 7b following the step of 
obtaining image data, the image data can be enhanced 104, This process can optionally 
be repeated 105, as desired. 

[086] Alternatively, as shown in Fig. 7c the derived quantitative and/or qualitative 

information can be compared to an image taken at T1 162, or any other time, if such 



24 



wo 2004/025541 



PCT/US2003/030004 



image is available. As shown in Fig. 7d the step of enhancing image data 104 can follow 
the step of obtaining image data 102. Again, the process can be repeated 105, as desired. 

[087] As previously described, some or all the processes illustrated in Figs. 7a-d 

can be repeated one or more times as desired. For example, locating a region of 
interest 100, obtaining image data 102, enhancing Image data 104, obtaining quantitative 
and/or qualitative data 120, administering candidate agent 150, comparing quantitative 
and/or qualitative information to a database 160, comparing quantitative and/or qualitative 
information to an image tal<en at a prior time, such as Ti, f 62, monitoring therapy 170, 
monitoring disease progress 172, predicting disease course 174 can be repeated one or 
more times if desired, 101, 103, 105, 121, 151, 161, 163, 171, 173, -/ZS respectively. 
Each of these steps can be repeated in one or more loops as shown in Fie. 7b, 176, 177, 
178, 179, 180, as desired or appropriate to enhance data collection. 

[088] Turning now to Fig. 8a, following the step of extracting image data from the 

RO1 102, the image can be transmitted 180. Transmission can be to another computer in 
the network or via the World Wide Web to another network. Following the step of 
transmitting the image 180, the image is converted to a pattern of normal and diseased 
tissue 190 Normal tissue includes the undamaged tissue located in the body part selected 
for study. Diseased tissue includes damaged tissue located in the body part selected for 
study. Diseased tissue can also include, or refer to, a lack of normal tissue in the body part 
selected for study. For example, damaged or missing cartilage would be considered 
diseased tissue. Once the Image is converted, it is analyzed 200. Fig. 8b illustrates the 
process shown in Fig. 8a with the additional step of enhancing image data 104. As will be 
appreciated by those of skill in the art, this process can be repeated 105 as desired. 

[089] As shown in Fig. 8c, the step of transmitting the image 180 illustrated in 

Fig. 8a is optional and need not be practiced under the invention. As will be appreciated 
by those of skill in the art, the image can also be analyzed prior to converting the image to 
a pattern of normal and diseased. Fig. 8d illustrates the process shown in Fig. 8c with the 
additional step of enhancing image data 104 which is optionally repeated 105, as desired. 

[090] As previously described, some or all the processes in Figs. 8a-d can be 

repeated one or more times as desired. For example, locating a region of Interest 100, 
obtaining image data 102, enhancing image data YO^.transmitting an image 180, 
converting the image to a pattern of normal and diseased 190, analyzing the converted 
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image 200, can be repeated one or more times if desired, 101 j 103, 105, 181, 191, 201 
respectively. 

[091] Fig. 9 shows two devices 900, 920 that are connected. Either the first or 

second device can develop a degeneration pattern from an image of a region of interest 
905. Similarly, either device can house a database for generating additional patterns or 
measurements 915. The first and second devices can communicate with each other in the 
process of analyzing an image, developing a degeneration pattern from a region of 
interest in the image, and creating a dataset of patterns or measurements or comparing 
the degeneration pattern to a database of patterns or measurements. However, all 
processes can be performed on one or more devices, as desired or necessary. 

[092] In this method the electronically generated, or digitized image or portions of 

the image can be electronically transferred from a transferring device to a receiving device 
located distant from the transferring device; receiving the transferred image at the distant 
location; converting the transferred image to a pattern of normal or diseased or abnormal 
tissue using one or more of the parameters specified in Table 1 , Table 2 or Table 3; and 
optionally transmitting the pattern to a site for analysis. As will be appreciated by those of 
skill in the art, the transferring device and receiving device can be located within the same 
room or the same building. The devices can be on a peer-to-peer network, or an intranet. 
Altematively, the devices can be separated by large distances and the information can be 
transferred by any suitable means of data transfer, including the World Wide Web and ftp 
protocols. 

[093] Alternatively, the method can comprise electronically transferring an 

electronically-generated image or portions of an image of a bone or a joint from a 
transferring device to a receiving device located distant from the transferring device; 
receiving the transferred image at the distant location; converting the transferred image to 
a degeneration pattern or a pattern of normal or diseased or abnormal tissue using one or 
more of the parameters specified in Table 1 , Table 2 or Table 3; and optionally 
transmitting the degeneration pattern or the pattern of normal or diseased or abnormal 
tissue to a site for analysis. 

[094] Another aspect of the invention is a kit for aiding in assessing the condition 

of a bone or a joint of a subject, which kit comprises a software program, which when 
installed and executed on a computer reads a degeneration pattern or a pattern of normal 



26 



wo 2004/025541 



PCT/US2003/030004 



or diseased or abnormal tissue derived using one or more of the parameters specified in 
Table 1 , Table 2 or Table 3 presented in a standard grapliics format and produces a 
computer readout. The kit can further include a database of measurements for use in 
calibrating or diagnosing the subject. One or more databases can be provided to enable 
the user to compare the results achieved for a specific subject against, for example, a 
wide variety of subjects, or a small subset of subjects having characteristics similar to the 
subject being studied. 

[095] A system is provided that includes (a) a device for electronically transferring 

a degeneration pattern or a pattern of normal, ^diseased or abnormal tissue for the bone or 
the joint to a receiving device located distant from the transferring device; (b) a device for 
receiving said pattern at the remote location; (c) a database accessible at the remote 
location for generating additional patterns or measurements for the bone or the joint of the 
human wherein the database includes a collection of subject patterns or data, for example 
of human bones or joints, which patterns or data are organized and can be accessed by 
reference to characteristics such as type of joint, gender, age, height, weight, bone size, 
type of movement, and distance of movement; (d) optionally a device for transmitting the 
correlated pattern back to the source of the degeneration pattern or pattern of normal, 
diseased or abnormal tissue. 

[096] Thus, the methods and systems described herein make use of collections of 

data sets of measurement values, for example measurements of bone structure and/or 
bone mineral density from x-ray images. Records can be formulated in spreadsheet-like 
format, for example including data attributes such as date of x-ray, patient age, sex, 
weight, current medications, geographic location, etc. The database formulations can 
further comprise the calculation of derived or calculated data points from one or more 
acquired data points, typically using the parameters listed in Tables 1 . 2 and 3 or 
combinations thereof. A variety of derived data points can be useful in providing 
information about individuals or groups during subsequent database manipulation, and are 
therefore typically included during database formulation. Derived data points include, but 
are not limited to the following: (1) maximum value, e.g. bone mineral density, determined 
for a selected region of bone or joint or in multiple samples from the same or different 
subjects; (2) minimum value, e.g. bone mineral density, determined for a selected region 
of bone or Joint or in multiple samples from the same or different subjects; (3) mean value, 
e.g. bone mineral density, determined for a selected region of bone or joint or in multiple 
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samples from the same or different subjects; (4) the number of measurements that are 
abnormally high or low, determined by comparing a given measurement data point with a 
selected value; and the like. Other derived data points include, but are not limited to the 
following: (1) maximum value of a selected bone structure parameter, determined for a 
selected region of bone or in multiple samples from the same or different subjects; (2) 
minimum value of a selected bone structure parameter, determined for a selected region 
of bone or in multiple samples from the same or different subjects; (3) mean value of a 
selected bone structure parameter, determined for a selected region of bone or in multiple 
samples from the same or different subjects; (4) the number of bone structure 
measurements that are abnormally high or low, determined by comparing a given 
measurement data point with a selected value; and the like. Other derived data points will 
be apparent to persons of ordinary skill in the art in light of the teachings of the present 
specification. The amount of available data and data derived from (or arrived at through 
analysis of) the original data provides an unprecedented amount of information that is very 
relevant to management of bone-related diseases such as osteoporosis. For example, by 
examining subjects over time, the efficacy of medications can be assessed. 

[097] Measurements and derived data points are collected and calculated, 

respectively, and can be associated with one or more data attributes to form a database. 
The amount of available data and data derived from (or arrived at through analysis of) the 
original data provide provides an unprecedented amount of information that is very 
relevant to management of musculoskeletal-related diseases such as osteoporosis or 
arthritis. For example, by examining subjects over time, the efficacy of medications can 
be assessed. 

[098] Data attributes can be automatically input with the electronic image and can 

include, for example, chronological information (e.g., DATE and TIME). Other such 
attributes can include, but are not limited to, the type of imager used, scanning 
information, digitizing information and the like. Alternatively, data attributes can be input 
by the subject and/or operator, for example subject identifiers, i.e. characteristics 
associated with a particular subject. These identifiers include but are not limited to the 
following: (1) a subject code (e.g., a numeric or alpha-numeric sequence); (2) 
demographic information such as race, gender and age; (3) physical characteristics such 
as weight, height and body mass index (BMI); (4) selected aspects of the subject's 
medical history (e.g., disease states or conditions, etc.); and (5) disease-associated 
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characteristics such as the type of bone disorder, if any; the type of medication used by 
the subject. In the practice of the present invention, each data point would typically be 
identified with the particular subject, as well as the demographic, etc. characteristic of that 
subject. 

[099] Other data attributes will be apparent to persons of ordinary skill in the art in 

light of the teachings of the present specification. (See, also, WO 02/30283, incorporated 
by reference in its entirety herein). 

[0100] Thus, data (e.g., bone structural information or bone mineral density 

information or articular information) is obtained from normal control subjects using the 
methods described herein. These databases are typically referred to as "reference 
databases" and can be used to aid analysis of any given subject's x-ray image, for 
example, by comparing the information obtained from the subject to the reference 
database. Generally, the information obtained from the normal control subjects will be 
averaged or otherwise statistically manipulated to provide a range of "normal" 
measurements. Suitable statistical manipulations and/or evaluations will be apparent to 
those of skill in the art in view of the teachings herein. The comparison of the subject's x- 
ray information to the reference database can be used to determine if the subject's bone 
information falls outside the normal range found in the reference database or is 
statistically significantly different from a normal control. 

[0101] Data obtained from images, as described above, can be manipulated, for 

example, using a variety of statistical analyses to produce useful information. Databases 
can be created or generated from the data collected for an individual, or for a group of 
individuals, over a defined period of time (e.g., days, months or years), from derived data, 
and from data attributes. 

[0102] For example, data can be aggregated, sorted, selected, sifted, clustered and 
segregated by means of the attributes associated with the data points. A number of data 
mining software exist which can be used to perform the desired manipulations. 

[0103] Relationships in various data can be directly queried and/or the data 

Analyzed by statistical methods to evaluate the information obtained from manipulating the 
database. 
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[0104] For example, a distribution curve can be established for a selected data set, 
and the mean, median and mode calculated therefor. Further, data spread characteristics, 
e.g., variabiiity, quartiles, and standard deviations can be calculated. 

[0105] The nature of the relationship between any variables of interest can be 
examined by calculating correlation coefficients. Useful methods for doing so include, but 
are not limited to: Pearson Product Moment Correlation and Spearman Rank Correlation. 
Analysis of variance permits testing of differences among sample groups to determine 
whether a selected variable has a discernible effect on the parameter being measured. 

[0106] Non-parametric tests can be used as a means of testing whether variations 
between empirical data and experimental expectancies are attributable to chance or to the 
variable or variables being examined. These include the Chi Square test, the Chi Square 
Goodness of Fit, the 2x2 Contingency Table, the Sign Test and the Phi Correlation 
Coefficient. Other tests include z-scores, T-scores or lifetime risk for arthritis, cartilage loss 
or osteoporotic fracture. 

[0107] There are numerous tools and analyses available in standard data mining 
software that can be applied to the analyses of the databases that can be created 
according to this invention. Such tools and analysis include, but are not limited to, cluster 
analysis, factor analysis, decision trees, neural networks, rule induction, data driven 
modeling, and data visualization. Some of the more complex methods of data mining 
techniques are used to discover relationships that are more empirical and data-driven, as 
opposed to theory driven, relationships, 

[0108] Statistical significance can be readily determined by those of skill in the art. 
The use of reference databases in the analysis of x-ray images facilitates that diagnosis, 
treatment and monitoring of bone conditions such as osteoporosis. 

[0109] For a general discussion of statistical methods applied to data analysis, see 

Applied Statistics for Science and Industry, by A. Romano, 1977, Allyn and Bacon, 
publisher. 

[0110] The data is preferably stored and manipulated using one or more computer 
programs or computer systems. These systems will typically have data storage capability 
(e.g., disk drives, tape storage, optical disks, etc.). Further, the computer systems can be 
networked or can be stand-alone systems. If networked, the computer system would be 
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able to transfer data to any device connected to the networked computer system for 
example a medical doctor or medical care facility using standard e-mail software, a central 
database using database query and update software (e.g., a data warehouse of data 
points, derived data, and data attributes obtained from a large number of subjects). 
Alternatively, a user could access from a doctor's office or medical facility, using any 
computer system with Internet access, to review historical data that can be useful for 
determining treatment. 

[01 1 1] If the networked computer system includes a World Wide Web application, 
the application includes the executable code required to generate database language 
statements, for example, SQL statements. Such executables typically include embedded 
SQL statements. The application further includes a configuration file that contains 
pointers and addresses to the various software entities that are located on the database 
server in addition to the different external and internal databases that are accessed in 
response to a user request. The configuration file also directs requests for database 
server resources to the appropriate hardware, as can be necessary if the database server 
is distributed over two or more different computers. 

[0112] As a person of skill in the art will appreciate, one or more of the parameters 
specified in Table 1, Table 2 and Table 3 can be used at an initial time point Ti to assess 
the severity of a bone disease such as osteoporosis or arthritis. The patient can then 
serve as their own control at a later time point Tz, when a subsequent measurement using 
one or more of the same parameters used at Ti is repeated. 

[01 13] A variety of data comparisons can be made that will facilitate drug discovery, 
efficacy, dosing, and comparisons. For example, one or more of the parameters specified 
In Table 1, Table 2 and Table 3 can be used to identify lead compounds during drug 
discovery. For example, different compounds can be tested in animal studies and the iead 
compounds with regard to highest therapeutic efficacy and lowest toxicity, e.g. to the bone 
or the cartilage, can be identified. Similar studies can be performed in human subjects, 
e.g. FDA phase I, II or III trials. .Alternatively, or in addition, one or more of the 
parameters specified in Table 1, Table 2 and Table 3 can be used to establish optimal 
dosing of a new compound. It will be appreciated also that one or more of the parameters 
specified in Table 1 .Table 2 and Table 3 can be used to compare a new drug against one 
or more established drugs or a placebo. 
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[01 14] The foregoing description of embodiments of the present invention has been 
provided for the purposes of illustration and description. It is not intended to be 
exhaustive or to limit the invention to the precise forms disclosed. Many modifications and 
variations will be apparent to the practitioner skilled in the art. The embodiments were 
chosen and described in order to best explain the principles of the invention and Its 
practical application, thereby enabling others skilled in the art to understand the invention 
and the various embodiments and with various modifications that are suited to the 
particular use contemplated. It is Intended that the scope of the invention be defined by 
the following claims and its equivalence. 
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What IS claimed is: 

1 . A method for analyzing at least one of bone mineral density, bone structure 
and surrounding tissue comprising: 

a. obtaining an image of a subject; 

b. locating a region of interest on the image; 

c. obtaining data from the region of interest; and 

d. deriving data selected from the group of qualitative and quantitative 
from the image data obtained at step c, 

2. A kit for aiding in the assessment of the condition of at least one of a bone 
and joint comprising: 

a software program which reads at least one of a degeneration pattern, a 
pattern of normal tissue, a pattern of abnormal tissue, and a pattern of diseased 
tissue. 

3. The kit of claim 2 further comprising a database of measurements for 
comparison to the at least one of degeneration pattern, pattern of normal tissue, pattern of 
abnormal tissue, and pattern of diseased tissue. 

4. The kit of claim 2 further comprising a subset of a database of 
measurements for comparison to the at least one of degeneration pattern, pattern of 
normal tissue, pattern of abnormal tissue, and pattern of diseased tissue. 

5. At least one of an automated and semi-automated method of using an 
imaging marker comprising: 

obtaining image data from a subject; 

obtaining data from the image data wherein the data obtained is at least one 
of quantitative and qualitative data; and 
administering an agent. 

6. The method of claim 5 wherein the one of the automated and semi- 
automated method is used for at least one of drug discovery, diagnosis, disease staging, 
disease monitoring, disease management, prognostication, therapy monitoring, drug 
efficacy monitoring, and disease prediction. 

7. At least one of a system for monitoring the efficacy of an agent and a system 
for drug discovery comprising: 

administering an agent to a subject; 
obtaining Image data; and 

obtaining data from the image data wherein the data obtained is at least one 
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of quantitative and qualitative data. 

8. At least one of a system for diagnosing a disease, determining disease 
staging, monitoring disease progression, managing a disease, disease prognostication, 
predicting a disease, monitoring therapy and randomizing a subject within a group of 
patients comprising: 

a. obtaining image data of a subject; 

b. obtaining data from the image data wherein the data obtained is at 
least one of quantitative and qualitative data; and 

c. comparing the at least one of quantitative and qualitative data in step 
b to at least one of: a database of at least one of quantitative and qualitative data 
obtained from a group of subjects; at least one of quantitative and qualitative data 
obtained from the subject; and at least one of a quantitative and qualitative data 
obtained from the subject at time Tn. 

9. The method of claims 1,5,7 and 8 further including the step of enhancing 
image data. 

10. The method of claims 1 , 5, 7 and 8 wherein the subject is a mammal. 

1 1 . The method of claims 1 , 5, 7 and 8 wherein the subject is a human. 

12. The method of claims 1 , 5, 7 and 8 wherein the subject Is a horse. 

13. The method of claims 1 , 5, 7 and 8 wherein the step of obtaining Image data 
Includes obtaining data from a measured parameter selected from the group consisting of: 
bone parameters, cartilage parameters, cartilage defect parameters, cartilage disease 
parameters, area parameters, and volume parameters. 

14. The method of claims 1 , 5, 7 and 8 wherein the step of obtaining image data 
includes extracting data from a measured parameter selected from the group consisting 
of: 

a presence or absence of bone marrow edema; 
a volume of bone marrow edema; 

a volume of bone marrow edema normalized by at least one of width, area, 
size, and volume; 

a presence or absence of osteophytes; 

a presence or absence of subchondral cysts; 

a presence or absence of subchondral sclerosis; 

a volume of osteophytes; a volume of subchondral cysts; 

a volume of subchondral sclerosis; 
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an area of bone marrow edema; an area of osteophytes; 

an area of subchondral cysts; 

an area of subchondral sclerosis; 

a depth of bone marrow edema; 

a depth of osteophytes; 

a depth of subchondral cysts; 

a depth of subchondral sclerosis; 

at least one of a volume, area, and depth of at least one of an osteophytes, 
subchondral cysts, subchondral sclerosis wherein the at least one of volume, area, 
and depth is normalized by at least one of width, area, size, volume a bone 
proximal to at least one of the osteophyte, subchondral cyst, or subchondral 
sclerosis; 

a presence or absence of meniscal tear; 
a presence or absence of cruciate ligament tear; 
a presence or absence of collateral ligament tear; 
a volume of menisci; 

a ratio of volume of normal to at least one of torn, damaged and 
degenerated meniscal tissue; 

a ratio of surface area of normal to at least one of torn, damaged and 
degenerated meniscal tissue; 

a ratio of surface area of normal to at least one of torn, damaged and 
degenerated meniscal tissue to total joint or cartilage surface area; 

a ratio of surface area of at least one of torn, damaged and degenerated 
meniscal tissue to a total surface area of at least one of joint and cartilage; 

a size ratio of opposing articular surfaces; 

a meniscal subluxation/dislocation in millimeters; 

an index combining different articular parameters ; 

a 3D surface contour information of subchondral bone; 

an actual or predicted knee flexion angle during gait cycle; 

a predicted knee rotation during gait cycle; 

a predicted knee displacement during gait cycle; 

a predicted load bearing line on cartilage surface during gait cycle and 
measurement of distance between load bearing line and at least one of cartilage 
defect and diseased cartilage; 
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a predicted load bearing area on cartilage surface during gait cycle and 
measurement of distance between load bearing area and at least one of cartilage 
defect and diseased cartilage; 

a predicted load bearing line on cartilage surface during standing or different 
degrees of knee flexion and extension and measurement of distance between load 
bearing line and at least one of cartilage defect and diseased cartilage; 

a predicted load bearing area on cartilage surface during standing or 
different degrees of knee flexion and extension and measurement of distance 
between load bearing area and at least one of cartilage defect and diseased 
cartilage; 

a ratio of load bearing area to area of at least one of cartilage defect and 
diseased cartilage; 

a percentage of load bearing area affected by cartilage disease; 

a location of cartilage defect within load bearing area; 

a load applied to cartilage defect, area of diseased cartilage; and 

a load applied to cartilage adjacent to at least one of cartilage defect and 
area of diseased cartilage. 

15. The method of claim 14 wherein the index combining different articular 
parameters includes: 

a presence or absence of cruciate or collateral ligament tear in the subject, 
a body mass index for the subject, 
a weight for the subject, or 
a height for the subject. 

16. The method of claims 1 , 5, 7 and 8 wherein the image data is obtained from 
a hip and the step of obtaining image data includes extracting data from a measured 
parameter selected from the group consisting of: 

microarchitecture parameters on structures parallel to stress lines; 
microarchitecture parameters on structures perpendicular to stress lines; 
geometry; 
shaft angle ; 
neck angle; 

diameter of femur neck; 
hip axis length; 

largest cross-section of femur head; 
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average thickness of cortical within at least one ROl; 
standard deviation of cortical thickness within at least one ROl; 
- maximum thickness of cortical within at least one ROl; 
minimum thickness of cortical within at least one ROl; and 
hip joint space width. 

17. The method of claims 1 , 5, 7 and 8 wherein the image data is obtained from 
a region of a spine and the step of obtaining image data includes extracting data from a 
measured parameter selected from the group consisting of: 

microarchitecture parameters on vertical structures; 
microarchitecture parameters on horizontal structures; 
geometry; 

superior endplate cortical thickness; 
inferior endplate cortical thickness; 
anterior vertebral wall cortical thickness; 
posterior vertebral wall cortical thickness; 
superior aspect of pedicle cortical thickness; 
inferior aspect of pedicle cortical thickness; 
vertebral height; 
vertebral diameter; 
pedicle thickness; 
, maximum vertebral height; 

minimum vertebral height; , 

average vertebral height; 

anterior vertebral height; 

medial vertebral height; 

posterior vertebral height; 

maximum inter-vertebral height; 

minimum, inter-vertebral height; and 

average inter-vertebral height. 

18. The method of claims 1 , 5, 7 and 8 wherein the image data is obtained from 
a region of a knee and the step of obtaining image data includes extracting data from a 
measured parameter selected from the group consisting of: 

average medial joint space width; 
minimum medial joint space width; 
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maximum medial joint space width; 

average lateral joint space width; 

minimum lateral joint space width; and 

maximum lateral joint space width. 
19. The method of claims 1 , 5, 7 and 8 wherein the step of obtaining Image data 
includes extracting bone parameters selected from the group consisting of: 

stainless steel equivalent thickness wherein the stainless steel equivalent 
thickness Is determined as the average gray value of the region of interest 
expressed as thickness of stainless steel with equivalent intensity; 

trabecular contrast wherein the trabecular contrast is determined as one of 
the trabecular equivalent thickness and marrow equivalent thickness; 

fractal dimension; 

Fourier spectral analysis wherein the Fourier spectral analysis is determined 
as one of a mean transform coefficient absolute value and a mean spatial first 
moment; 

predominant orientation of spatial energy spectrum; 
at least one of trabecular area and total area; 
trabecular perimeter; 
trabecular distance transform; 
marrow distance transform; 

trabecular distance transform regional maxima values; 
marrow distance transform regional maxima values; 
star volume; 

trabecular bone pattern factor ; 

connected skeleton count or trees (T); 

node count (N); 

segment count (S); 

node-to-node segment count (NN); 

node-to-free-end segment count (NF); 

node-to-node segment length (NNL) 

node-to-free-end segment length (NFL); 

free-end-to-free-end segment length (FFL); 

node-to-node total struts length (NN.TSL); 

free-end-to-free-ends total struts length( FF.TSL); 
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total struts length (TSL); 
FF.TSL/TSL; 
NNTSL/TSL; 
loop count (Lo); 
loop area; 

mean distance transform values for each connected skeleton; 
mean distance transform values for each segment (Tb.Th ); 
mean distance transform values for each node-to-node segment (Tb.Th.NN); 
mean distance transform values for each node-to-free-end segment 
(Tb.Th.NF); 

orientation of each segment; 
angle of each segment; 
angle between segments; 

length-thickness ratios (NNL/Tb.Th.NN ) and (NFU Tb.Th.NF); and 
interconnectivity index (ICI) ICI = (N * NN)/ ( T * (NF + 1) ); 

20. The method of claim 1 2 wherein total bone parameter factor is (P1 - P2) / 
(A1 - A2 ) 

further wherein P1 and A1 are the perimeter length and trabecular bone area 
before dilation and P2 and A2 corresponding values after a single pixel dilation, 
measure of connectivity 

21 . The method of claims 1,5,7 and 8 wherein the step of obtaining image data 
includes extracting at least one of cartilage parameters, cartilage defect parameters, and 
diseased cartilage parameters, wherein the data extracted is selected from the group 
consisting of: 

total cartilage volume; 
focal cartilage volume; 

a cartilage thickness distribution or thickness map; 
mean cartilage thickness over substantially total surface; 
mean cartilage thickness in focal area; 
median cartilage thickness; 
maximum cartilage thickness; 
minimum cartilage thickness; 
3D cartilage surface information; 
cartilage curvature analysis; 
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volume of cartilage defect/diseased cartilage; 
depth of cartilage defect/diseased cartilage; 
area of cartilage defect/diseased cartilage; 

at least one of 2D and 3D location of cartilage defect/diseased cartilage in 
the articular surface; 

at least one of 2D and 3D location of cartilage defect/diseased cartilage in 
relationship to weight-bearing area; 

a ratio of at least two of diameter of cartilage defect, diameter of diseased 
cartilage, and thickness of surrounding normal cartilage; 

a ratio of at least two of depth of cartilage defect, depth of diseased 
cartilage and thicl<ness of surrounding normal cartilage; 

a ratio of at least two of volume of cartilage defect, volume of diseased 
cartilage and thickness of surrounding normal cartilage; 

a ratio of at least two of surface area of cartilage defect, surface area of 
diseased cartilage and total joint surface area; and 

a ratio of at least two of volume of cartilage defect, volume of diseased 
cartilage and total cartilage volume, 

22. The method of claims 1,5,7 and 8 wherein the steps are performed 
automatically. 

23. The method of claims 1,5,7 and 8 wherein the steps are performed semi- 
automatically. 

24. The method of claims 1,5,7 and 8 wherein at least one of the steps of the 
method is performed on a first computer. 

25. The method of claims 1,5,7 and 8 wherein at least one of the steps of the 
method is performed on a first computer and at least one of the steps of the method is 
performed on a second computer. 

26. The method of claim 25 wherein the first computer and the second computer 
are connected. 

27. The method of claim 26 wherein the first computer and the second computer 
are connected by one of a peer to peer network, direct link, intranet, and internet. 

28. The method of claim 1 wherein the step of locating a region of interest is 
repeated. 

29. The method of claims 1, 5, 7 and 8 wherein the step of obtaining image data 
from a region of interest is repeated. 
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30. The method of claims 1, 5, 7 and 8 wherein at least one of the image and 
image data is converted to a 2D pattern. 

31 . The method of claim 30 wherein the 2D pattern is evaluated. 

32. The method of claim 30 wherein the 2D pattern is converted to a 3D pattern. 

33. The method of claim 31 wherein the 2D pattern is converted to a 3D pattern. 

34. The method of claim 30 wherein the 2D pattern is converted to a 4D pattern. 

35. The method of claim 31 wherein the 2D pattern is converted to a 4D pattern. 

36. The method of claim 33 wherein the 3D pattern is converted to a 4D pattern. 

37. The method of claims 1, 5, 7 and 8 wherein at least one of the image and 
image data is converted to a 3D pattern. 

38. The method of claim 37 wherein the 3D pattern is evaluated. 

39. The method of claim 37 wherein the 3D pattern is converted to a 2D pattern. 

40. The method of claim 39 wherein the 3D pattern is converted to a 2D pattern. 

41 . The method of claim 37 wherein the 3D pattern is converted to a 4D pattern. 

42. The method of claim 38 wherein the 3D pattern is converted to a 4D pattern. 

43. The method of claim 39 wherein the 2D pattern is converted to a 4D pattern. 

44. The method of claims 1,5,7 and 8 wherein at least one of the image and 
image data is converted to a 4D pattern. 

45. The method of claim 44 wherein the 4D pattern is evaluated. 

46. The method of claim 1 further comprising the step of administering a 
candidate agent. 

47. The method of claim 46 wherein the candidate agent is at least one agent 
selected from the group consisting of: substance administered to a subject, substance 
ingested by a subject, molecules, pharmaceuticals, biopharmaceuticals, 
agropharmaceuticals. 

48. The method of claims 1, 5, 7 and 8 further comprising the step of comparing 
at least one of the image and image data to a database. 

49. The method of claims 1,5,7 and 8 further comprising the step of comparing 
the at least one of the image and image data to a subset of a database. 

50. The method of claims 1, 5, 7 and 8 further comprising the step of comparing 
at least one of quantitative data and qualitative data to an image taken at T1 . 

51 . The method of claims 1,5,7 and 8 further comprising the step of comparing 
at least one of quantitative data and qualitative data to an image taken prior to the image 
under analysis. 



41 



wo 2004/025541 



PCT/US2003/030004 



52. The method of claims 1,5,7 and 8 further comprising the step of comparing 
at least one of quantitative data and qualitative data to an image taken at Tn. 

53. The method of claim 1 further comprising the step of transmitting at least one 
of the image, data extracted from the region of interest. 

54. The method of claims 1 , 5, 7 and 8 wherein at least one of the image and 
image data is converted to at least one of a pattern of normal, a pattern of diseased, and a 
pattern of normal and diseased. 

55. The method of claims 1 , 5, 7 and 8 wherein at least one of obtaining image 
and image data includes measuring at least one of microarchitecture and 
macroanatomical structures. 

56. The method of claim 55 further comprising measuring the average density. 

57. The method of claim 56 wherein the average density measurement includes 
a calibrated density of the region of interest. 

58. The method of claim 55 further comprising measuring microanatomical 
structures on at least one of dental, spine, hip, knee and bone core x-rays. 

59. The method of claim 58 further comprising measuring at least one of: 
calibrated density of extracted structures; 

calibrated density of background; 
average intensity of extracted structures; 

average intensity of background area wherein the background area includes 
non-extracted structures; 

structural contract wherein structural contrast is an average intensity of 
extracted structures divided by an average intensity of a background; 

calibrated structural contrast wherein calibrated structural contrast is a 
calibrated density of extracted structures divided by a calibrated density of a 
background; 

total area of extracted structures; 

total area of a region of interest; 

an area of extracted structures normalized by a total area of a region of 
interest; 

a boundary length of an extracted area normalized by a total are of a region 

of interest; 

a number of structures normalized by an area of a region of interest; 
a trabecular bone pattern factor; 
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a measurement of concavity and convexity of structures; 
a star volume of extracted structures; 
a star volume of bacl^ground; and 

a number of loops normalized by an area of a region of interest. 

60. The method of claim 55 further comprising measuring a distance transform 
of extracted structures. 

61 . The method of claim 60 wherein the measurement on the distance transform 
of extracted structures further comprises one or more of: 

an average regional maximum thlcl^ness; 
a standard deviation of regional maximum thickness; 
a largest value of regional maximum thickness; and 
a median regional maximum thickness. 
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